TeV Gamma Rays from Geminga and the Origin of the GeV Positron Excess 
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The Geminga pulsar has long been one of the most intriguing MeV-GeV gamma-ray point sources. 
We examine the implications of the recent Milagro detection of extended, mufti- TeV gamma-ray 
emission from Geminga, finding that this reveals the existence of an ancient, powerful cosmic-ray 
accelerator that can plausibly account for the multi-GeV positron excess that has evaded explana- 
tion. We explore a number of testable predictions for gamma-ray and electron/positron experiments 
(up to ~ 100 TeV) that can confirm the first "direct" detection of a cosmic-ray source. 

PACS numbers: 95.85.Ry, 98.70.Rz, 98.70.-f 



Introduction. — Geminga holds a place of distinction 
among gamma-ray sources, being the first pulsar to be 
discovered through gamma rays, with a history of obser- 
vations through a variety of techniques While one 
of the brightest MeV-GeV gamma-ray point sources in 
the sky, there was no certain evidence of high-energy ac- 
tivity beyond the immediate neighborhood of the pulsar 
or its x-ray pulsar wind nebula (PWN) until the recent 
detection by Milagro of gamma rays at ~ 20 TeV from 
a region of ~ 3° around the pulsar [3, 0] . This detection 
places Geminga among the growing class of TeV PWNe 
(e.g., and is important for understanding aged pul- 

sars and their winds. An immediate consequence is the 
existence of a population of high-energy particles. 

The relative proximity of Geminga raises an interesting 
possibility, namely that these high-energy particles, most 
likely electrons and positrons, may be at the root of the 
explanation of the "positron excess" , the observed @, 0, 
[1| overabundance of multi-GeV positrons as compared 
to theoretical expectations (see Fig. [1]) . Severe energy 
losses of high-energy positrons require a local source of 
some kind [lfj, such as Geminga lllor even dark matter 
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through its annihilation products 

Here, we connect the Milagro TeV gamma-ray "halo" 
to electrons and positrons with energies up to at least 
100 TeV, expected to be accelerated in PWNe (e.g., 
13l . 14 1; for a review see 15|), and present several pre- 
dictions. Principally, while Geminga is apparently young 
enough to still produce high-energy particles, it is old 
enough that multi-GeV electrons and positrons from its 
more active past could have made it to Earth. The ex- 
tended gamma-ray emission is strong evidence for pro- 
duction, acceleration, and escape, suggesting an explana- 
tion of the positron excess. Moreover, this single nearest 
high-energy astrophysical source can reasonably account 
for the e~ + e + spectrum as measured by Fermi [16| and 
HESS 
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18l ] with an extension to energies beyond sev- 
eral TeV, where no signal might be expected otherwise. 

The Gamma-ray Source Next Door. — The observation 



of high-energy gamma rays from an astrophysical source 
implies the presence of higher-energy particles, typically 
or protons, that gave rise to them. One striking ele- 
ment of the observation of ~ 20 TeV gamma rays (with a 
significance of 4.9 a in the PSF-smoothed map [2j, 6.3 a 
for an extended source [3j) from Geminga by Milagro is 
the extent of the emission, 9 ~ 3° Q , which corresponds 
to a physical size of sq ~ 10 pc (9g/3° )(rc/200pc), 
where rQ is the distance to Geminga. Since the angu- 
lar resolution of Milagro is better than a degree and the 
characteristic age of the pulsar, to ~ 3 x 10 5 yr [l~9j ]. 
seemingly excludes a typical TeV supernova remnant, we 
shall consider an extended PWN with emission from a 
much larger region than seen in x-rays [2(J|2lj]. We will 
draw guidance from the TeV-PWN HESS J1825-137 0], 
which, while only a tenth the age of Geminga, would 
appear tens of degrees wide if placed at tq ~ 200 pc. 

We first examine whether the gamma rays can be ex- 
plained through inverse-Compton (IC) up-scattering of 
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FIG. 1: The cosmic-ray pos itron fraction. Shown are data 
compiled from Refs. @, S HI [34j , and scenarios based on the 
secondary model of Ref. [9 1 (shaded) and a plausible Geminga 
contribution (solid, dashed, and dotted lines) dependent upon 
distance and energetics (see text for details). 
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cosmic microwave background (CMB) photons by e ± . 
Note that the pulsar age exceeds the IC cooling time on 
CMB photons of the > 100 TeV e ± needed to produce 
>20 TeV gamma rays, t ic ~ 10 4 (100 TcV/£7 e ) yr in the 
Thomson limit. Including synchrotron losses further de- 
creases t coo i, implying fresh e ± production. To account 
for the Milagro measurement of 6.9 ± 1.6 x 10~ 15 TeV -1 
cm~ 2 s~ 4 at 20 TeV (see Fig. [2]), we consider a generic 
parent spectrum of the form dN/dj oc ^~ Q e ~ 7 / 7max , 
with 7 = E/(m e c 2 ). Lacking more detailed observa- 
tions, we choose E min = 1 GeV, 2? max = 200 TeV and 
a = 2 (typical to shock acceleration and as inferred in 
the Vela X PWN [1]). The resulting IC spectrum is 
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/f dl^J 
dj J dE b — n b (E b ) cr KN (-f, E b , £7 7 ), (1) 



where the Klcin-Nishina cross section, okn, is given by 
Ref. [22|]. The dominant scattering background, n b (E b ), 
in the Milagro energy range is the CMB, allowing us to 
construct the minimal spectrum shown in Fig. [5] (at lower 
energies, other contributions become relevant). When we 
normalize the IC spectrum to the Milagro TeV gamma- 
ray luminosity of £ 7 ,TcV ~ 10 32 erg s" 1 (for tq ~ 200 pc), 
at least £ e ± ~ 10 45 erg of is required. 

We can relate this phenomenological spectrum to 
the pulsar. The Goldreich-Julian flux [13] is Nqj ~ 
B tt 2 R 3 /ec. For Geminga, R ~ 10 km, B ~ 1.6 x 10 12 G, 
and n ~ 26.5 s" 1 0, which yield A> GJ ~ 10 32 s _1 . 
For our E~ 2 spectrum, (E e ±) w 15 GeV. The total 
power is then £ e ± w N e ± (E e ±) rj MN GJ (E e ±) w 
2 x 10 30 M erg s -1 , where M = N e ± /iV G j is the pair mul- 
tiplicity. For a pure electron flow (M ~ 1), the power 
beyond 10 TeV is only - 5 x 10 29 erg s" 1 < £ 7 , To v, 
thus requiring pair production resulting in M. ^ 100. 
Including synchrotron losses comparable to IC increases 
the required A4. This can only be avoided by assuming 
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FIG. 2: Minimal inverse-Compton gamma-ray spectrum of 
the extended emission from Geminga (shaded) and the Mi- 
lagro measurement at 20 TeV (left axis). Also, the energy 
distribution (dotted line) of the associated e (right axis). 



a spectrum much-harder than E~ 2 or E m i n > 100 GeV 
(both inconsistent with, e.g., Vela X @, H^]). 

While we have been trying to explain the observed sig- 
nal alone, there may well be more lower surface bright- 
ness emission at larger angles, and also at lower energies. 
Two considerations make this likely. First, the pulsar's 
e^- output was probably much stronger in the past when 
its spin-down power was higher. Secondly, no evidence of 
a large-scale radio or x-ray nebula exists, and a substan- 
tial fraction of the e may be escaping, so that the above 
multiplicity is only a lower limit. Ai ~ 10 (as inferred 
for younger TeV PWNe [23] or pulsar models [24Jj ) does 
not exceed the spin-down power of ~ 10 34 5 erg s _1 [lj|, 
so that a large pair conversion fraction is possible. 

The Origin of the Positron Excess. — The confirmed 
presence of a nearby, ancient source of high-energy elec- 
trons and positrons immediately suggests an explanation 
for the positron excess. If so, then we would essentially 
be living "within" the extended halo of the source, seeing 
that were accelerated long ago when the pulsar was 
stronger. The density of particles at a given time and 
place, n(r, t, 7), is governed by the diffusion equation, 
which (in spherically symmetric geometry) is 
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where Q is the source term. In our range of interest, 
the energy loss rate is well approximated by £(7) = lol 2 ■ 
Considering IC losses on the CMB (energy density ~ 
0.3 eVcm -3 ) and synchrotron losses due to a ~ 3 /iG 
magnetic field (~ 0.2 eV cm -3 ) yields 4 - 5x 10~ 20 s~ 4 . 

Calculations that assess the positron flux from dark 
matter annihilation typically assume steady-state con- 
ditions (i.e., ignoring the time-dependent term). How- 
ever, our scenario is manifestly dynamic, with both the 
source luminosity and distance potentially changing with 
time. The analytic solution given by Atoyan, Aharonian, 
and Volk 25j allows for this to be conveniently handled. 
Their diffusion coefficient is 2? (7) = X>o(l + 7/7* ) 5 , with 
7, ~ 6 x 10 3 [25)]. In the limit of a single burst from a 
point source, Q("f) — dN / dry 8(r) S(t — to), with 

n(r t ^ - W. [r/rS^fe-^n)? 
n\r, 1, 7; — /2 M n x _ A2 _ ra ■ \o) 



(l-^7) 2 



where the spectrum is cut off at j c = 1/ (£q t) . The energy 
loss rate and age then set the maximal energy of particles 
that reach us today, with a diffusion radius of 



r d (t, 7) ~ 2 (2%) t [1 - (1 - 7/7c) 1 -"]/[(l - 5)1 he]) 
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For ( G ~3x 10 5 yr, 2? n 



' yr, L>o ~ 4 x 10 7 cm s~\ and 5 = 0.4 
(intermediate between 6=1/3 and 1/2 [26j), the diffu- 
sion radius is r d ~ 150, 175, 250 pc for £7 = 2, 10, 50 GeV 
particles, respectively, so multi-GeV particles are now ar- 
riving. We caution against extrapolating to small radius, 
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FIG. 3: Cosmic-ray e" + e+ data from Fermi Iff. HESS [U 
GjJ, and Refs. [HEH and AMS e"-only data 01 with Galac- 
tic e~ model modified from Ref. [§[ (solid), contributions from 
Geminga (lower lines), and total (upper lines). 



since the diffusion solution may transition [2 



like 



form (as does HESS J1825-137 



to a wind- 
) near the 



source, although the data allow no firm conclusions. 

For a continuously emitting source such as Geminga, 
the injection rate can be parametrized as dN/dj oc 
£e± (t)l~ ae ~ 7 / 7max , with C e ± the luminosity. The lo- 
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cal particle density is 71-0(7) = J dtn(rG,t,j). Assum- 
ing braking via magnetic dipole radiation, the spin-down 
luminosity evolves as oc (1 + t/to)~ 2 [28|, with a pulsar- 
dependent timescale, to, and C e ± (t) = (£gAg) [1 + (to — 
t)/t ]- 2 /f G dt'[l + (t G - t')/t ]-\ For t a ~ 3 x 10 4 yr, 
the present spin-down power, ~ 10 34 ' 5 erg s _1 , corre- 
sponds to an upper limit on the total output of 
~ 5 x 10 48 erg (larger for smaller tn [25|). Geminga's 
transverse velocity is ~ 200 km s _1 [29(. A similar radial 
velocity would result in a ~ 100 pc displacement in ta- 
in Fig. [3j we display the local flux of + e + , 
Jq — (c/iir) uq, from our benchmark model of a = 2, 
within a reasonable range of parameters. These have dis- 
tances varying (from birth — > present) as r G — 150 — > 
250 pc, 220 pc, 250 -> 200 pc; e ± energy budgets of 
£ G = 1 , 2 , 3 x 10 48 erg; and 5 = 0.4, 0.5 , 0.6, respec- 
tively (lower dotted, solid, dashed lines). The energy in 
e ± estimated for several younger TeV PWN are at least 
as large as these (e.g., [J, [5J, |23|). Since the bulk of the 
energy is released in this early spin-down phase, the ini- 
tial location is the most important. Adding to these the 
primary e~ spectrum of Moskalenko and Strong with 
the normalization decreased by 35% and an added ex- 
ponential cutoff at 2 TeV (in order to not exceed HESS 
data), yields the total e~ + e + flux (upper lines). 

The spectral feature at ~ 1 TeV naturally results from 
a combination of energy losses and pulsar age and dis- 
tance (see Fig. 4 of Ref. [ljj for comparison) . The multi- 



TeV extension (beyond the last HESS point) is due to the 
continuous injection of particles, as evidenced by the Mi- 
lagro observations today. Combining these with the ex- 
pectations for the secondary e ± fluxes (see also 
we compare our positron fraction to measurements in 
Fig- HI (note that solar modulation may account for dis- 
agreements between data below ~ 10 GeV l3ll]). 

It is thus plausible that Geminga is the long-sought [32[ 
local source of electrons and positrons, influencing the 
spectra measured by Fermi [16j (down to tens of GeV) 
and HESS 
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in the TeV, although we emphasize 
that certain parameters and the underlying Galactic pri- 
mary spectrum remain uncertain. The PAMELA [33j 



and AMS |34| experiments can measure the e~ and e + 
spectra separately to isolate this component (since the 
e~ spectrum from Geminga should be identical to the 
e+). 

Conclusions. — The discovery of high-energy gamma 
rays from an extended region around Geminga by Mila- 
gro reveals the presence of > 100 TeV e ± , as observed 
indirectly within the x-ray PWN 2^, 21 1 . A considerable 
amount of data should become available as new experi- 
ments examine the surrounding area. This will help in 
developing more detailed models that account for both 
time and spatial evolution in the e spectra, directly 
coupled to cosmic-ray propagation [35j. One need is a 



better-determined distance, the most recent quoted be- 
ing tq ~ 250lg2° P c [29'J ■ We briefly discuss implications 
for several categories of experiments. 

Fermi: While the observed features of Geminga will 
depend upon details such as whether the source is roughly 
spherical or preferentially oriented, we would generally 
expect the source to become "larger" with decreasing 
energy, reflecting the decrease in IC cooling time with 
energy. Our inspection of the point-source subtracted 
sky map from EGRET [36j indicates emission in the 
GeV range of a size comparable to the Milagro source. 
Fermi [37j should be able to more effectively separate the 
bright pulsed signal to study diffuse emission. 

Te V gamma rays: Obtaining a detailed spectrum and 
morphology of the source in the TeV regime will be vi- 
tal for further interpretation of the nature of the parti- 
cles present. Already, VERITAS [H] has placed rather- 
tight upper limits on a point source at the location of 
Geminga [i^]. Further study of the expected extended 
source is needed to better estimate the total energet- 
ics. In HESS J1825-137, the surface brightness was seen 
to drop off as ~ 1/0, inconsistent with pure diffusion 
and suggestive of convection, and the gamma-ray spec- 
trum was measured to soften with increasing distance 
from its pulsar Q. We expect similar behavior from 
Geminga if the same mechanisms are at work, the lat- 
ter of which would be a distinct signature of e ± cool- 
ing Q. Also, studying the extended TeV emission from 
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an old, radio-quiet neutron star should have implications 
for some heretofore unidentified TeV sources. 

Electrons and Positrons: Due to the spin down of the 
pulsar, it is possible that the Geminga source was much 
brighter in the past and dominated the TeV sky. It is 
from this time that multi-GeV e may still be reaching us 
today. If Geminga does account for a substantial fraction 
of the total e~ + e + spectrum, a mild anisotropy may be 
present [13, 41 1. Since the distance to Geminga does 
not gre atly exceed the scale for field fluctuations of ~ 
100 pc 26], with detailed multi- wavelength studies, local 
diffusion parameters might be determined (which may 
differ from global values estimated across the Galaxy). 
Additionally, as Geminga remains a source of ~ 100 TeV 
e ± , it may result in a > 10 TeV lepton flux at Earth. 

Neutrinos: The gamma rays might be produced via 
the decay of neutral pions produced in hadronic scatter- 
ing [42| if a nucleonic wind carries away the spin-down 
energy (as proposed for Vela X; but see Ref. [23j]). 
However, this requires ~ 10 47-48 erg of protons to be 
present and confined for > 10 5 yr (likely confining as 
well), disfavoring this scenario. There would also be a 
neutrino flux from charged pions. This would be simi- 



lar to MGRO J20 19+37 [44J, so that IceCube [45( would 
expect > 1 (0.2) neutrino- induced muon per year with en- 
ergy > 1 (10) TeV (see 46]). An improved measurement 
of the source extent is needed to estimate the atmospheric 
background, which can exceed this rate for a radius much 
larger than a degree [3] . Even if no neutrinos are found, 
IceCube will be able to provide valuable constraints 47 1 . 
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